Suppression effect on explosive percolations 



O 
(N 

o 

CD 

Q 



o 
o 

O 
• I— I 

Oh: 



> 
o 

00 
O 



X 




Y.S. Cho and B. Kahng 
Department of Physics and Astronomy, Seoul National University, Seoul 151-747, Korea 

(Dated: December 30, 2011) 

Percolation transitions (PTs) of networks, leading to the formation of a macroscopic cluster, are 
conventionally considered to be continuous transitions. However, a modified version of the classical 
random graph model was introduced in which the growth of clusters was suppressed, and a FT 
occurs explosively at a delayed transition point. Whether the explosive PT is indeed discontinuous 
or continuous becomes controversial. Here, we show that the behavior of the explosive PT depends 
on detailed dynamic rules. Thus, when dynamic rules are designed to suppress the growth of all 
clusters, the discontinuity of the order parameter tends to a finite value as the system size increases, 
indicating that the explosive PT could be discontinuous. 

PACS numbers: 02.50.Ey,64.60.ah,89.75.Hc 



Percolation transition (PT), i. e., the transition from 
a disconnected state to a connected one, has been re- 
garded as a fundamental model of phase transitions in 
noncquilibrium systems [l| . The concept of PT has been 
extended to the formation of macroscopic clusters in net- 
work science. A pioneering model of PT in network sci- 
ence is the classical random graph model introduced by 
Erdos and Renyi (ER) 0] in which a system composed of 
a fixed number of vertices N evolves as edges are added. 
At each evolution step, an edge is added between two 
vertices, that are randomly selected from among uncon- 
nected pairs of vertices. In this model, a quantity, called 
time, is defined as the number of edges added to the sys- 
tem per node. The ER model has been modified by fol- 
lowing the so-called Achlioptas process [i]. The Achliop- 
tas process essentially identifies the dynamics that pre- 
vent the creation of a given target pattern by choosing 
one edge from a given number of randomly selected po- 
tential edges. In the modified ER model for the PT, 
the target pattern is a giant cluster. Thus, an edge to 
be added to the system should be selected such that the 
growth of clusters can be systematically suppressed. The 
principle to take this selection rule is hereafter referred 
to as the suppression principle (SP). 

The dynamic rule originally designed by Achlioptas et 
al. 0] is as follows: in the case of two randomly selected 
edge candidates, the one actually added to the system 
is the one minimizing the product or the sum of the 
sizes of the clusters that are connected by each potential 
edge. The ER models modified according to the prod- 
uct rule and the sum rule are denoted as ERPR and 
ERSR models, respectively. In these models, the giant- 
cluster size increases drastically at the critical point, and 
therefore, the percolation transition is called explosive 
percolation. The introduction of an explosive PT model 
has triggered intensive research on discontinuous PTs in 
non-equilibrium systems Ii-12|. Many models have fol- 
lowed the ERPR and ERSR models, and they display 
similar transition patterns. Although the explosive PT 
was regarded as discontinuous in the original paper 



FIG. 1. (Color online) Classification of types of edge can- 
didate pairs, (i) Both candidates ei and 62 are intercluster 
edges, (ii) ei is an intracluster edge and the other candidate 
62 is an intercluster edge, (iii) Both ei and 62 are intracluster 
edges. 



recently it has been argued that the transition is continu- 
ous in the thermodynamic limit 



13l-ll7l|. Thus, the issue 



of whether the explosive PT is indeed discontinuous or 
continuous remains controversial. In this Letter, we de- 
scribe the microscopic investigation of the dynamic rules 
of several explosive percolation models and the surpris- 
ing finding that these rules do not satisfy the SP. Thus, 
it is rather natural that the PTs of those models are con- 
tinuous. However, some other variants of the Achlioptas 
model satisfying the SP exhibit the pattern of discontin- 
uous PTs within the range of our numerical simulations. 
Thus, we may state that satisfying the SP is essential for 
discontinuous PTs in the evolution of complex networks. 

To explain the dynamic rule, we classify the types of 
edge candidate pairs as follows: 

(i) Both edge candidates ei and 62 arc intercluster 
edges. Clusters of sizes sJ2 ^-iid ^il 
nected by the edge ei, and clusters of sizes and 

are connected by the edge 62. In the follow- 
ing discussion, we will use the following notation: 
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(ii) ei is an intracluster edge in a cluster of size sj"\ 
and 62 is an intercluster edge between two clusters 
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(iii) Both ei and 62 are intracluster edges in either the 
same cluster or two different clusters. 

These three types of edge candidate pairs are depicted in 
Fig. [TJ On the basis of this classification, we formulated 
several dynamic rules to determine which edge should be 
added to the system. 

Although the original model Q seems to follow the 
basic idea of the Achlioptas process, the dynamic rule 
has to be more carefully examined. As time approaches 
the percolation threshold, the mean cluster size increases, 
and one or both potential edges have a greater possibility 
of being intracluster as shown in Fig. [H^a). Thus, we 
have to clarify how to formulate the dynamic rule when 
intracluster edges are selected as potential edges. Here, 
we formulate dynamic rules for the cases (i)-(iii), and we 
check whether each rule does follow the SP. 

First, we introduce three different variants of the 
ERFR model; these are specified in Table L In model A, 
when one edge is an intracluster edge ((ii) and (iii)), the 
product is the square of the size of that cluster, while for 
case (i), it is the product of the sizes of the two clusters 
connected by one intercluster edge. This rule, however, 
can fail to follow the Achlioptas SP. For example, when 
sj"^ = 5, Sja = 3, and 325^ = 7 for case (ii), the edge 
62 is selected in the ERPR-A model, and then, the size 
of the created cluster is 10. On the other hand, if edge 
ei is selected, then none of the clusters would increase in 
size. Therefore, model A does not follow the SP. As the 
transition point is approached, intracluster edges in (ii) 
and (iii) can be selected more frequently [see Fig. [21Ja)]. 
Actually, the behavior of the fraction of the occurrences 
of (ii) or (iii) is similar to that of G{t). Thus the failure 
of the SP can be more frequent. 

In the ERPR-B model, when potential edges of the 
type (ii) are selected, the intracluster edge is definitely 
selected, so that the cluster size does not increase. When 
the two potential edges are both intracluster edges (type 
(iii)), one of them is randomly selected. By this rule, 
time is advanced by one unit 1/iV for the types (ii) and 
(iii). 

Model C is a simplified version of models A and B. In 
this model; the dynamics proceed via only intercluster 
connections. Thus, two potential edges are both inter- 
cluster edges (case (i)). Model C may be regarded to be 
nearly the same as model B because the clusters do not 
grow when the intracluster edge is selected in (ii) and 
(iii). However, the difference between them is that for 
types (ii) and (iii), time is advanced in model B but not 
in model C. 

For all models A, B, and C, the product rule has an 
intrinsic drawback that the Achlioptas SP is unfulfilled. 
Let us consider a simple example of two intercluster con- 
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FIG. 2. (Color online) (a) The fractions of type (ii) and type 
(iii) potential edges as a function of time t for different system 
sizes: iV/10'' = 32, 64, 128, and 256 from the top (bottom) 
in the small-f (large-f) region. As increases, the fraction 
increases dramatically. Inset: For a given N — 1.024 x lO'^, 
the fractions of the type (i) (dotted line) and the types (ii) 
and (iii) (solid lines) as a function of t. They are compared 
with the giant-cluster size Gjv(t) (dashed line). The fraction 
of the types (ii) or (iii) behaves similarly to Gjv(t), indicating 
that counting for the effect by taking the intracluster edges 
becomes important as t approaches tc- (b) The fraction of 
the occurrences in which the sum of the sizes of one pair of 
clusters becomes larger than that of the size of the other pair, 
even though the product of the sizes of the former is smaller 
than that of the sizes of the latter. The dotted line represents 
Gjv(t). Inset: Solid lines represent the failure ratio on an 
enlarged scale for different system sizes N/W^ = 32,64, 128, 
and 256. For a larger system, the curve lies on the upper 
position in the small-t region. The dashed line represents 
Gjv(t). Numerical data for (a) and (b) are obtained from the 
ERPR-B model. 



nections, in which s^l — 2 



- 7, 41 = 4, and s^^ 



Then, P, 



(i) 



14 and P: 



(i) 



'Ih ~ ' ' ''2a — '^26 

16, and thus, edge Ci is added 
to the system. However, the resulting cluster size is 9 in 
the case, which is larger than the resulting size 8 when 
62 is added. In other words, even though the product of 
one pair of cluster sizes is smaller than that of the other 
pair, its sum can be larger. Thus, the Achlioptas SP is 
inherently unfulfilled. Investigations 0, [13, 13, 16| have 
shown that the cluster size distribution displays a hump 
shape in the region of large cluster sizes, and the hump 
size increases up to the point where explosive cluster ag- 
gregations start. Owing to the inherent drawback, such 
a case is likely to occur frequently, as shown in Fig.[2Jb). 
Thus, the PTs under the product rule are continuous re- 
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TABLE I. List of the dynamic rules under the product rule [the sum rule] . The second column lists the type of potential edges. 
The third and fourth columns list the conditions for the cases in which both ei and 62 in Fig. 1 are selected. The last column 
shows the case when either ei or 62 is selected randomly. 



Model 



Type 



ei 



£2 



Either ei or £2 randomly 



Model A 



Type (i) 
Type (ii) 
Type (iii) 



(i) 
2 



■5f ) < S«i 

(ii) 



> P^" [S[' 



(i) 



> S. 



{i)l 



p(i) 



3(i) 
2 

D{ii) 



(ii)l 



j(i)i 



,(ii)). 



D(ii) 



[24" 



(ii) 



?(ii)l 



unconditional 



Model B 



Type (i) 
Type (ii) 
Type (iii) 



unconditional 



^(i)l 



' > p^'> [55" > s. 



(i) tqW 



?{i)l 



P^' 



(i) 



3(i) rQ(i) 



?(i)l 



unconditional 



Model C 



Type (i) 



P« < P^w [S[^> < S^' 



(i) rQ{i) 



j{i)l 



P« > P^W [51" > 5^' 



(i) rQ(i) 



?{i)l 



P,^" = P,« [5« = 



?(i)l 



gardless of the model type. 

Next, wc introduce similar models under the sum rule; 
these are also listed in Table I. The drawback inherent to 
the product rule is removed in the ERSR model. How- 
ever, for case (ii), the Achlioptas SP cannot be fulfilled 
in model A, but it is always fulfilled in models B and C. 
Thus, in the case of the sum rule, the models B and C 
are regarded to be the ones following the Achlioptas SP. 
We state that the PTs for models B and C are possible 
candidates for the discontinuous PT. 

Even though it is a challenging task to determine the 
transition types of the explosive PTs with numerical sim- 
ulation data, the numerical approach is the only one pos- 
sible, since there is no analytic solution that takes into ac- 
count all the aforementioned cases. Extensive numerical 
simulations are carried out up to a system size TV = 10^" 
with a configuration average of about The ob- 

tained numerical data may be sufficient for understand- 
ing why this controversy has arisen, but a higher configu- 
ration average may be required for determining the type 
of PT, particularly when the system size is large. 

We measure GN{t) as a fimction of time t for different 
system sizes in the range iV = 2° x 10* — 2^° x 10'' at 
every N = 2 x 10^ step. For given N and 2N, we find a 
point of intersection of the two curves G'jv(i) and G2jv(i)- 
The time and G components of such a point are denoted 
as t^{2N) and G2Ar(ix(2A^)), respectively. We compose a 
set of {tx{N), GN{tx{N))} for the simulated system sizes. 
To evaluate the discontinuity of the PT, we propose the 
following criteria: 

(a) The value (tx(-^), GAr(tx)) remains finite as N — > 
00. The time tx(oo) is regarded as the transition 
point tc in the thermodynamic limit. 

(f3) The tangent of the curve GN(t) with respect to t at 
t-^{N) diverges as N increases. 

Figs. [3Ua)-(c) show the behaviors of GAr(tx) as a func- 
tion of N for models A, B, and C, respectively, under 



the product rule. The insets of each figure show GAr(t) 
versus t. We can see in (a) and (c) that for the model 
A and C under the product rule, GAr(tx) decreases with 
increasing N in the whole considered range, suggesting 
that GN{tc) — !• in the limit TV — > cx). In the case of 
model B, even though the data look flat up to ~ 10*, 
they decay in the large- A'^ region in the same manner as 
for model C. Thus, the decay behavior can be consid- 
ered to stem from the intrinsic drawback of the product 
rule for the case (i). In our previous study [lo| . we per- 
formed numerical simulations for the model B under the 
product rule up to a system size N = 10*. In that case, 
however, the decay behavior was not noticed and the PT 
was regarded as a discontinuous transition. With the 
simulation data obtained in this sutdy for larger system 
sizes, we conclude that the three models based on the 
product rule show continuous PTs. 

Figs. [3l[d)-(f) are the plots of GNitx) versus N for 
models A, B, and C under the sum rule. It can easily be 
seen that for model A, the GN{tx) values decrease with 
increasing N, suggesting that GM{tc) is zero in the limit 
A — !■ 00. This indicates that the rule of doubling the 
cluster size in the ERSR-A model violates the SP, and 
leads to a continuous transition. However, for models B 
and C, the data of the GN{tx) values look relatively fiat 
asymptotically within the large- A^ limit, even though the 
data points beyond A = 10^ have large error bars owing 
to a smaller number of configuration averages. Fig. [2Jg) 
shows the tangent to the curve GAr(t) at the crossing 
point as a function of A^. The data show that the tangent 
increases according to a power law ~ A^'^'^. Thus, we may 
conclude that the ERSR-B and ERSR-C models fulfill the 
SP and seem to show discontinuous PTs, within the range 
of our numerical data. However, we cannot rule out the 
possibility that the value of Gjv(tx(A)) decreases when 
the system sizes are larger than those simulated in this 
work. 

Recently, the authors [l3| introduced a new type of 
Achlioptas percolation model and argued that the ex- 
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FIG. 3. (Color online) G]v(tx) versus iV for the (a) ERPR- 
A, (b) ERPR-B, (c) ERPR-C, (d) ERSR-A, (e) ERSR-B, 
and (f) ERSR-C models. The slopes of each guideline is 
(a) -0.06, (b) -0.05, (c) -0.05, (d) -0.02, (e) 0, and (f) 
0. Thus, Giv(tx) of ERSR-B and ERSR-C converge to fi- 
nite values in the thermodynamic limit. The error bars rep- 
resent the deviation of the cross points. Each data set is 
obtained by taking an average over about 10^^ /N config- 
urations. Gjv(tx) of models B and C of both the ERPR 
and ERSR overlap in the large-A?^ region. The insets of 
(a)-(f) show the behaviors of Gjv(t) for different system 
sizes iV/lO" = 32, 64, 128, 256, 512, 1024, 2048, 4096, 8192, and 
16384. (g) Plot of the slope of G]v(t) at versus A'^ for the 
ERSR-B (□) and the ERSR-C (o) models. dGN{t^)/dt in- 
creases according to a power law ~ N°'^. 



plosive PT is actually continuous. This model, called 
the CDGM model following the initials of the authors, 
is defined as follows: Firstly, a pair of clusters, Cia and 
Ci6, are randomly selected and the smaller cluster (say 
Cia) is selected. Secondly, another pair of clusters C2a 
and 6*26 are randomly picked, and the smaller one (say 
C2a) is selected. Thirdly, two random nodes from each of 
the chosen clusters (Cia and C2a) are selected and con- 
nected. There are four possible combinations of the con- 
nection. However, when either of the two clusters from 
the first set is identical to cither of the two clusters from 
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FIG. 4. (Color online) (a) Schematic illustration of the se- 
lection rule on the intracluster edge for the CDGM-B model. 
In the original CDGM model (CDGM-A), when two pairs of 
edges ei and e-i are selected from the cluster sets Gi and G2 
of sizes sia < sii, and S2a < S2b, an additional edge connects 
the two nodes in the clusters C\a and C^a (dotted line). In 
our modified model (CDMA-B) the two nodes in the clus- 
ters Gib and C2h are connected instead (solid line), because 
by choosing them, the cluster sizes in the system do not in- 
crease. Panels (b)-(d) show GAr(fx) versus A'^ for CDGM-A, 
CDGM-B, and CDGM-C models, respectively, and the slopes 
are —0.03, and respectively. GAr(tx) of the models B and C 
overlap in the large A region. It may be reasonable to expect 
that the modified CDGM model (CDGM-B and C) shows a 
discontinuous PT. 



the second set (for example, C\b and C2h in Fig. |4l[a)), 
the CDGM model can fail to follow the Achlioptas SP. 
According to the original rule of the CDGM model, the 
two smaller clusters C\a and C^a are connected, which 
creates a larger cluster size whereas the connection be- 
tween two nodes inside the cluster C\h = C^b does not 
increase any cluster size in the system. Therefore, the 
original CDGM model fails to follow the Achlioptas SP 
by not taking into account the natural choice of that 
intra-cluster edge. Again, as time approaches tc, the se- 
lection of an intracluster edge becomes more frequent. 
Thus, the selection of intercluster edges, which is against 
the SP, can change the PT into a continuous transition. 
We confirm this result by performing the following nu- 
merical simulations. We plot GN{tx) versus TV, and find 
Gw(tx) - for the CDGM model (Fig. Hb)). We 

then modify the original model as follows: when one or 
two intracluster edges are present among the four edge 
candidates, one of these intracluster candidates is con- 
nected (model B) randomly. In addition, we consider a 
model C in which the four edge candidates are only al- 
lowed to be intercluster edges. Figs.lDJc) and (d) suggests 
that GAr(<x) approaches a constant value in the large- A'^ 
region for models B and C, respectively. Moreover, we 
confirm that the slope of GN{t) at diverges as A^ in- 
creases in a power law manner. Based on these numeri- 
cal results, it reveals that the PT could be discontinuous 
even in the CDGM model, when properly modified. 

In summary, we have examined the dynamic rule of 
the Achlioptas model in the perspective of the supprcs- 
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sion principle (SP). We found that when the dynamic 
rule does follow the SP, the numerically calculated order 
parameter seems to show the behavior of a discontinuous 
transition. Otherwise, the PT seems to be continuous. 
The original Achlioptas model and the CDGM model be- 
long to the latter case. An analytic study of the modified 
Achlioptas models by taking into account the suppression 
effect is therefore needed. 

This study was supported by an NRF grant awarded 
through the Acceleration Research Program (Grant No. 
2010-0015066) and the NAP of KRCF (BK) and the 
Seoul Science Foundation (YSC). 
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